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Abstract. Observations have revealed strong evidence for powerful jets in the Low/Hard states of black hole 
candidate X-ray binaries. Correlations, both temporal and spectral, between the radio - infrared and X-ray bands 
suggest that jet synchrotron as well as inverse Compton emission could also be significantly contributing at higher 
frequencies. We show here that, for reasonable assumptions about the jet physical parameters, the broadband 
spectrum from radio through X-rays can be almost entirely fit by synchrotron emission. We explore a relatively 
simple model for a relativistic, adiabatically expanding jet combined with a truncated thermal disk conjoined by 
an ADAF, in the context of the recently discovered black hole binary XTE J1118+480. In particular, the X-ray 
power-law emission can be explained as optically thin synchrotron emission from a shock acceleration region in 
the innermost part of the jet, with a cutoff determined by cooling losses. For synchrotron cooling-limited particle 
acceleration, the spectral cutoff is a function only of dimensionless plasma parameters and thus should be around a 
"canonical" value for sources with similar plasma properties. It is therefore possible that non-thermal jet emission 
is important for XTE J1118+480 and possibly other X-ray binaries in the Low/Hard state. 
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1. Introduction 

Observations are providing increasing evidence that black 
hole candidate (BHC) X-ray binaries (XRBs) produce 
powerful collimated outflows when in the Low/Hard X- 
ray state (LHS). This state is characterized by a nonther- 
mal power-law in the X-ray band and little, if any, ther- 
mal disk contribution (e.g., Nowak 1995| ; Poutancn 1998| ). 
At least three systems in the LHS (Cyg X-l, IE 1740.7- 
2942 and GRS 1758-258) have directly resolved radio jets 
on scales from AU to parsecs. However, XRB jets also 
reveal themselves in the broadband LHS spectra with a 
flat-to-inverted radio synchrotron spectrum, analogous to 
the sig nature emission of jets in compact radio cores of 
AGN ([Blandiord fc Konigl 1979j ; |Hjellming fc Johnston 



198S ; Falcke Sz Biermann 1999). Furthermore, new data 



suggest a continuation of this radio synchrotron emission 
to much higher frequencies. The spatial, spectral and tem- 
poral evidence for powerful jets fro m XRB BHC s in the 
LHS are compiled and discussed in Fender (2001 ). 



We know that jets play a significant role in the emis- 
sion of Active Galactic Nuclei (AGN), even dominating 
the spectrum from radio through TeV 7-rays in the case 
of BL Lacs, with emission from optically thin synchrotron 



up to even 100 keV and higher (e.g. Pian ct al. 1998). 
By analogy, if the flat, optically thick synchrotron spec- 
trum in XRBs, commonly attributed to j ets, indeed ex - 
tends into the NIR and optical regimes ( Fender 2001 ), 
one would expect a corresponding optically thin power- 
law from shock acceleration at even higher frequencies. 
This is irrespective of whether the jet emission is boosted 
or not, since optically thick and optically thin emission 
would have similar Doppler factors. Shock acceleration is 
likely to be present in XRBs, given that optically thin 
power-law spectra are observed during their radio out- 
bursts (e.g. 



Fender & Kuulkers 2001, and refs. therein) 



Still, the majority of current models for the broadband 
(X-ray) spectra of BHC XRBs focus only on the contri- 
bution of thermal disk plu s coronal invers e Compton (IC) 



emission (for a review, see Poutanen 1998 ) and ignore any 
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As an example for the possible importance of jet emis- 
s ion, we use the recent ly discovered XRB XTE J1118+480 
( Remillard et al. 200C| ) , which has been observed in the ra- 



dio through X-rays (see Hynes et al. 2000, hereafter H00 



Fender et al. 2001, hereafter F01, and references therein), 
and is also at high enough Galactic latitude to allow the 
first ever EUV detections of an X-ray transient (H00). 



The system is likely a BHC in the LHS (H00;R.cvnivtsev 



Sunyae |v, fc Borozdin 2000|; |Wood et al. 2000[). Although 
jets were not directly resolved with MERLIN to a limit 
of < 65(<i/kpc) AU (at 5 GHz; F01), its radio emission 
shows the fiat characteristic jet spectrum. 

2. Basic Model and Estimates 

We consider that an accretion disk is responsible for the 
extreme-ultraviolet (EUV) data and contributes in the op- 
tical (H00; |G arcia et al. 2000]). However, the power-law 



to be of order gj = lO^-lO" 3 ( |Falcke fc Biermann 1999| ). 
However, q } is essentially a free parameter, with q } <C 1, 
and we define q h ~2 = <ft/10 . While the jet formation 
itself is very difficult to model, the physics of calculating 
most of the jet emission is relatively straightforward be- 
cause the fiat-to-inverted spectrum stems from the part 
of the jet where it is basically undergoing free expansion. 
Here we build on the jet emission model outlined in Falcke 
fc Markoff (2000] ), and references therein 



seen ir the X-ray spectrum indicates that an additional 
component must be present. So far, this has been pre- 
sumed to arise from the Comptonization of some 'seed' 
photons by a hot, thermal corona (Poutanen 1998). 



At the inner edge of the hot accretion flow, plasma is 
ejected out from symmetric nozzles, where it becomes su- 
personic. The jets then accelerate along the axes through 
their pressure gradients up to bulk Lorentz factors 7j ~ 
2 — 3, and expand sideways with their initial proper sound 
speed 7s/3 s c — 0.4c for a hot electron/proton plasma. This 
implies low relativistic Mach numbers around M. ~ 5. 
The velocity field, density and magnetic field gradient then 
come naturally from the Euler equation (see, e.g., Falcke 



1996). The dependencies of the magnetic field B and den- 



sity n on distance are then similar to, but slightly stronger 
than, the canonical r _1 and r~ 2 dependencies for conical 



LHS is that a standard thin, optically thick disk (Shakura 
fc Sunyaev 1973) exists only down to some transition 
radius r tr ~ 10 2 — 10 3 r s , (r s = 2GMbh /c 2 ), where the 



jets, respectively (Blandford & Konigl 1979; Hjellming & 



Onj; commonly invoked physical explanation tor the Johnston 1988|; [Falcke fc Biermann 1995|) 



flow becomes ho t and non-radiative (e.g., |Liu et al. 1999] ; 
Esin et al. 2001 ). We do not, however, include an entire 



In this way, the basic physical properties governing the 
emission at each point in the jet are fixed after specify- 
ing the jet power, and the initial conditions at the noz- 
zle. We make the simplification of assuming a maximal 
jet, which follows from the Bernoulli equation (Falcke fc 



m existence explo ring this issue (see above, and review 
in Poutanen 1998 ). Instead, we include a representative 
thermal spectrum for the inner edge of the cool disk both 
in the direct emission, in our calculations of cooling rates 
in the jet and as seed photons for IC. 

With this approach the temperature, T d , an d lumi- 
nosity, L d , for this inner edge can be roughly determined 
by fitting a black body to the EUV data. As discussed 
in H00, uncertainty in the local absorption leads to large 
variations in the possible EUV flux for XTE J1118+480. 
Here, we take the highest absorption value presented, 
Ah = 1-15 • 10 20 cm -2 , because it provides a solid up- 
per limit to any thermal disk contribution (see Fig. 1). 

Taking the distance and BH mass to be 1.8 kpc and 
6M Q , respectively (McClintock et al. 2001a), the fit to the 
EUV data gives T d ~ 1.5 • 10 5 K and L d ~ 5 ■ 10 35 erg 
s _1 . For an annulus with scale width ~ r tr , we find as 
an order of magnitude estimate r tr ~ 900r s , in agreement 
with current models (e.g., Liu et al. 1999] ). The luminosity 
of a standard accretion disk is dominated by the inner edge 
and its radiati ve efficiency is q\ = jr s /r tI (Frank, King, 
fc Raine 1992 , Eq. 5.20), yielding a rough estimate for 
the accretion rate of M d ~ q^L d /c 2 ~ 3 ■ 10" 8 M Q yr" 1 . 
Hence, in the following we use a reference value of M d = 
m_7.5 3 x 10~ 8 Af Q yr _1 . Within r tr we consider a hot, 
ADAF-likc flow which does not significantly contribute to 
the spectrum. 

For accreting black holes it has been argued that the jet 
power is of order Qj Ct ~ qjM d c 2 with an efficiency inferred 



disk mjodcl here, since there are several models already Biermann 1995|) when the internal energy, here dominated 



by the magnetic field, is equal to the bulk kinetic energy 
of particles. This is consistent with a magnetic launching 
mechanism. The plasma is assumed to originate in the 
hot accretion flow and therefore contains equal numbers 
of protons and electrons, with hot electrons at a temper- 
ature appro aching T c = T e , io • 10 10 K in various ADAF 
models (e.g., Manmoto 2000). The electron Lorentz factor 
of the peak will be at 7 C ~ 4 • T e io. 

In AGN jets, the high frequency, optically-thin power- 
laws are taken to be the result of synchrotron emission 
from particles being sh ock-accelerated along the jet (e.g., 
Marscher fc Gear 1985 ) . In such a case the crucial param- 
eter for the high energy emission is the location z acc of the 
first particle acceleration region in the jet. 

Near the shock region in each jet, B 2 (r)/8ir ~ 
0.25gjMdC 2 /(c7j/3j7rr 2 ), yielding a reference value of B ~ 
2 • 10 6 G • \J gj j _2rh_7.5/7j/3j(r/10r s ) _1 for the parameters 
discussed above. Similarly, the particle density is given by 
n ~ 10 14 cm- 3 • ( % , _ 2 m_ 7 .5/7j/3j)(^/10r s )- 2 . 

Once the plasma, assumed to be injected at the base 
of the jet with a Maxwellian distribution, reaches the 
shock region, the standard diffusive shock acceleration 
process redistributes the particles into a power-law, start- 
ing roughly at 7 c , m in ~ 4T C! io. The initially injected dis- 
tribution has an index p~ 1.5 — 2 (resulting fro m a rela- 
tivistic shock, see, e.g., Heavens fc Drury 198Sf ) steepen- 



1 Because the plasma is mildly relativistic at the shock re- 
gion, the temperature-dependent adiabatic index will be closer 
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ing by ~ 1 due to increased cooling above a break energy 
E\y « (T r /3 coo i) _1 , where r r is the residence time of the 
plasma in the acceleration region, and /3 coo i are the con- 
stants giving the energy dependence of the cooling, which 
for synchrotron and Compton losses go as E ss —/3 coo iE 2 . 
This steepened spectrum has the index of p ~ 2.5 — 3, 
typically found in the optically thin synchrotron emission 
of both AGN and X-ray binaries. For the case of XTE 
J1118+480, the unbroken X-ray power-law (see Fig. 1), 
implies p ~ 2.6 for E > E^. The acceleration ceases when 
the particles reach the energy -E^max = r Ye,max'mcC 2 where 
the cooling/loss rates equal that of acceleration. These 
rates are dependent both on the energy of the particle, as 



well as the local physical parameters. 



in the IR/optical regime at ~ 10 14 5 Hz. For the parame- 
ters discussed here this would be at a distance of roughly 
z acc ~ 50r s , with a jet radius of r ~ 10r s , and the mag- 
netic field and particle densities derived above. The total 
synchrotron luminosity integrated up to the highest en- 
ergies, is a few times 10 36 erg s _1 , which is ~ 10% of the 
total jet power (Qj ~ <?j, -2^-7.5 2 • 10 37 erg s _1 ). 

Some fraction of X-rays created close to the nucleus 
by the jet will of course either directly impinge on, or be 
scattered by hot electrons into, the cold disk, resulting 
in a reflection component. However, since this is very de- 
pendent on the disk geometry, and because the feature 
is very weak or absent for XTE J1118+480 (|McClintock 



The shock acceleration rate is given as 



et al. 2001b| ), we do not calculate this process here. 



t: 



eB 



m c c£,j c 



(1) 3. Results of Numerical Modeling 



where w s h is the shock speed in the plasma frame. The 
parameter £ < c/3 e /u s h ( Jokipii 1987 ) is the ratio between 
the diffusive scattering mean free path and the gyroradius 
of the particle, and has a lower limit at £ = 1. We ac- 
count for energy losses at the shock via adiabatic losses, 
particle escape, IC and synchrotron. In our case the latter 
dominates and we have 



"syn 



4 r2 U b 

= t;Ot7oPc 

6 m e c 



(2) 



where <tt is the Thomson cross-section and Ub is the en- 
ergy density of the magnetic field. The maximum energy 
is then found by solving t^ c = t^, yielding 7 , max ~ 

10 8 (£-B) -0 ' 5 (^r)- If we define as a reference value £ = 
^100, the maximum synchrotron frequency is 



= 0.29^ r ~ 1.2 • 10 



20 1- 1 f Ush 
42 V c 



Hz 



(3) 



where v c ~ J^iUax 



(eB) I (m e c) is the critical synchrotron 
frequency. This maximum corresponds approximately to 
the rollover of the power-law cutoff, and for u a \ ~ j3 s c, we 
find a cutoff of ~ 80 keV. This cutoff is not dependent 
on the magnetic field, the jet power, or the shock location 
as long as we are in the synchrotron cooling dominated 
regime. Because we would expect XRBs to have similar 
shock structures, once £ is roughly fixed observationally — 
thus determining the scattering between magnetic irregu- 
larities in the diffusive shock process — we should get sim- 
ilar cutoffs for different sources and accretion rates. 

The location of the initial shock acceleration re- 
gion is set by the frequency where the flat, highly self- 
absorbed synchrotron spectrum turns over into the opti- 
cally thin power-law produced at the shock. From back- 
extrapolating the X-ray power-law, we can see that this 
maximum self-absorption frequency has to be somewhere 



to 4/3 rather than the non-relativistic value of 5/3. This yields 
a spectral index of 1.5 5s p < 2, where p = 2 is for non- 
relativistic shock-acceleration. 



To obtain a more detailed broadband spectrum we have 
used the full numerical calculations for a jet model as de- 
scribed in Falckc fc Markoff (2000] ), with the addition of a 
particle acceleration region. The model takes into account 
the relativistic Doppler shifts, adiabatic losses, electron 
acceleration and loss timescales at the acceleration region 
as described above, and integrates the synchrotron and IC 
emission along the jet. The resulting fit is shown in Figure 
1 with the parameters given in the caption. 

The relevant free parameters are the inclination angle 
9i, the shock distance z acc , the scattering ratio £, and the 
jet power parameterized by q y While the electron temper- 
ature T c (or 7c,min) and the fraction of thermal particles 
accelerated can also be adjusted, they are not indepen- 
dent of (J - decreasing these parameters will decrease the 
radiative efficiency an d increase qj . Sim ilarly, a more real- 
istic disk model (e.g., Esin et aL 2001) could yield a lower 
rtr and M, also increasing qj. Because we require charge 
balance between protons and electrons for the jet plasma, 
and a low T c in line with ADAF models, the particles are 
in sub-equipartition with the magnetic field. 

The flat spectrum from radio to optical is due to 
the optically thick synchrotron emission from the jet at 
z > z acc , while the X-ray power-law is the optically thin 
synchrotron emission dominated by the shock acceleration 
region (at z ~ z acc ). For the conditions in the jet, E\> lies 
very close to the peak of the thermal electron distribu- 
tion, and we only see the steepened power-law. Because 
of the large ratio between magnetic field and photon den- 
sities, IC emission does not play a significant role for the 
parameters chosen here (see Fig. 1). Synchrotron and IC 
emission from the pre-shock region (z < z acc ) could in 
principle show up in the soft and very hard X-rays respec- 
tively as a function of T e , but the EUV points place an 
effective upper limit T e <,2- 10 10 K, and so for this source 
these features are rather weak. 

The jet length may be constrained by the turnover 
seen in the ~ 2 — 15 GHz range, which gives a length of 
~ 4- 10 13 cm, i.e. ~ 1.5 mas (at 1 .8 kpc), in th e GHz range 
with a i^ -0 - 9 scaling of the size ( Falckc 1996| ). 
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injecting (mildly) relativistic electrons into the jet and 
allowing high electron energies because of low IC cool- 
ing. With a more dominant disk contribution, e.g., in the 
High/Soft state, the jet spectrum would also change and 
possibly disappear due to increased cooling. Interestingly, 
in the case of a low disk contribution, synchrotron- 
dominated cooling provides a natural "canonical" cutoff 
around 100 keV, which only depends on the dimensionless 
plasma parameters £ and u s h/c- Small variations in the 
spectral index of the synchrotron power-law could explain 



the observed tim e lags in some sources (Kotov, Churazov, 



fc Gilfanov 200l|) . The jet model for XTE J1118+480 may 



5 10 
log 10 v (GHz) 

Fig. 1. Fit to data from H00, with the exception of the 
VLA radio points (Hjellming et al., private comm.), the 
Ryle Telescope and a non-simultaneous JCMT point (see 
F01, and refs. therein). Included is a schematic indicating 
the different emission components. Parameters are Mbh = 
6M , r nozz = 1.5r s , T d = 1.7 • 10 5 K, L d = 5.25 • 10 35 
ergs s" 1 , m = 0.1, L } = 2.6 • 10 36 erg s -1 , r tr = 747r s , 
q it -2 = 1.5, £2 = 1, z acc = 45r s , 9 t = 32°. The fraction 
of thermal particles accelerated at the shock was fixed at 
50%. The dot-dashed line shows the cutoff for the case of 

6 = 0.1 



The optically thin synchrotron power-law component 
in the X-rays depends on the existence of a diffusive shock 
acceleration region. The current limit on y G comes from 
synchrotron losses instead of IC, due to the weak disk 
emission, which results in a low external photon density. 
In this case, the location of the X-ray cutoff is determined 
by £ and is fit observationally, keeping in mind that £ is 
typically under a few 10 2 (e.g., Jokipii 1987). A "canoni- 
cal" cutoff at ~ 100 keV (e.g., Poutancn 199S) corresponds 
to £ ~ 100, and u max oc We also plot both £ = 10 
and £ = 100 in Figure 1, in order to illustrate how the 
cutoff could theoretically venture into the realm of future 
high-energy missions like INTEGRAL with energies a 
few MeV. For X TE J1118+480, such a hi gh cutoff may 
be necessary, as VlcClintock et al. (2001b ) see no cutoff 
below ~ 160 keV. 



4. Conclusion 

From our modeling we conclude that in the LHS of 
XTE J1118+480 a significant jet contribution from radio 
through IR is present, and with reasonable assumptions 
and a small efficiency could extend even up to the hard 
X-rays. We cannot exclude that XTE J1118+480, with 
its rather low luminosity, is a special source. However, our 
model is scalable and thus likely applicable to other XRBs 
in the LHS, when a flat radio spectrum is present. 

An important element of the model is that the in- 
ner AD AF- like flow is very hot and radiatively inefficient, 



et al., 2001, 



2001, MNRAS, in press (astro- 



thcrefore provide an interesting new perspective for the 
modeling of XRBs in general. 
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